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ABSTRACT 

We present the discovery of the highest velocity C iv broad absorption line to date 
in the z=2.47 quasar SDSS J023011.28+005913.6, hereafter J0230. In comparing 
the public DR7 and DR9 spectra of J0230, we discovered an emerging broad 
absorption trough outflowing at ~60,000 kms“^, which we refer to as trough A. 
In pursuing follow up observations of trough A, we discovered a second emergent 
Civ broad absorption trough outflowing at ^40,000 kms“^, namely trough B. In 
total, we collected seven spectral epochs of J0230 that demonstrate emergent and 
rapidly (^10 days in the rest-frame) varying broad absorption. We investigate 
two possible scenarios that could cause these rapid changes: bulk motion and 
ionization variability. Given our multi-epoch data, we were able to rule out some 
simple models of bulk motion, but have proposed two more realistic models to 
explain the variability of both troughs. Trough A is likely an augmented ‘crossing 
disk’ scenario with the absorber moving at 10,000 < v (kms“^) < 18,000. Trough 
B can be explained by a flow-tube feature travelling across the emitting region at 
8, 000 < V (kms“^) < 56,000. If ionization variability is the cause for the changes 
observed, trough A’s absorber has rig > 724 cm“^ and is at Tequai ^ 2.00 kpc, or 
is at r < 2.00 kpc with no constraint on the density; trough B’s absorber either 
has He ^ 1540 cm“^ and is at Tequai ^ 1-37 kpc, or is at r < 1.37 kpc with no 
constraint on the density. 

Key words: quasars: general - quasars: absorption lines - quasars: emission 
lines - quasars: super massive black holes - quasars: individual: SDSS J0230+0059 
accretion discs 


1 INTRODUCTION 

At least 23% of quasars exhibit blueshifted Broad Absorp¬ 
tion Line (BAL) troughs at ultraviolet wavelengths (see 
discussions in Rogerson et al. 2011 and Allen et al. 2011), 
and the fraction increases if narrower (500-2000 km s“^) 
‘mini-BAL’ troughs are included (see Rodriguez Hidalgo 
et al. 2011 for a full discussion on mini-BAL quasars).^ 

* E-mail: rogerson@yorku.ca 

^ BAL quasars are, historically, defined as quasars that exhibit 
blueshifted absorption due to the Civ doublet at AA 1548.203, 
1550.770 A that is at least 2000 km s”*^ wide and can extend 
from 3000 kms“^ to 25000 kms“^, where 0 kms“^ is at the 
systemic redshift of the quasar (Weymann et al. 1991) and 
positive velocities indicate motion toward the observer. Modi- 


The disk-wind model of luminous Active Galactic Nuclei 
(AGN) characterizes BAL features as a result of mate¬ 
rial lifted off the accretion disk surrounding the central 
supermassive black hole (SMBH) and accelerated by ra¬ 
diation line driving to high outflow velocities that we ob¬ 
serve as blueshifted absorption (e.g. Murray et al. 1995, 
Ostriker et al. 2010). Whatever their origin, quasar out¬ 
flows provide insight into the physical and chemical prop¬ 
erties of the central engine of quasars, and may also rep¬ 


lications to this definition have been proposed, e.g. Hall et al. 
(2002) and Trump et al. (2006). Herein, we consider absorption 
at any velocity to be a candidate BAL trough, and we report 
widths of confirmed troughs so that others may classify the 
troughs as they see fit. 
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resent a mechanism by which SMBHs provide feedback to 
their host galaxy (e.g., Moe et al. 2009, Arav et al. 2013, 
Leighly et al. 2014, Chamberlain et al. 2015). 

Variability in the strength (i.e., the depth, width, or 
outflow velocity profile) of BALs is a well documented 
phenomenon (e.g., Gibson et al. 2008, Hall et al. 2011, 
Filiz Ak et al. 2013, He et al. 2015). Specifically, there 
have been recent studies documenting the disappearance 
of BAL troughs (e.g., Filiz Ak et al. 2012) as well as emer¬ 
gence in quasars that were not classified as having BALs 
previously (e.g., Rodriguez Hidalgo et ah, in preparation, 
Hamann et al. 2008, Leighly et al. 2009). 

The cause of BAL-trough variability is still largely 
debated in the literature, however, it is likely either due 
to transverse motion of absorbing clouds across our line 
of sight (e.g., Hall et al. 2011, Muzahid et al. 2015), or 
due to changes in the ionization of the absorbing gas 
(e.g., Hamann et al. 2008, Filiz Ak et al. 2013, Wang 
et al. 2015). Ultimately, it may be a complex mixture of 
these two scenarios. Full characterization of BAL variabil¬ 
ity events (either emergence, disappearance, or variability 
in general) would significantly increase our understanding 
of both the physics of the quasar’s central engine and the 
interaction of the quasar with its host galaxy. 

In this work we present the discovery of 
the highest velocity outflow discovered to date 
(~60,000 kms“^) at ultraviolet wavelengths, in the 
quasar SDSS J023011.28-b005913.6, hereafter J0230 
(Schneider et al. 2007). The previous highest-velocity 
absorption identified at ultraviolet wavelengths in a 
BAL quasar was at 56,000 kms“^ in PG 2302-1-029 
(Jannuzi et al. 1996), with the next highest being at 
50,000 kms“^ in PG 0935-1-417 (Rodriguez Hidalgo et al. 
2011).^ Outflows at these extremely high velocities have 
been previously observed in X-rays (e.g., Ghartas et al. 
2002, Pounds et al. 2003, but see Zoghbi et al. 2015) and 
might pose problems for theoretical acceleration models. 

We adopt a redshift oi z = 2.473 ± 0.001 for J0230 
based on visual inspection of the Lya, Gill], and Mgll 
emission lines and the onset of the Lya forest. Our red- 
shift is identical within the errors to the value of a = 
2.4721 ±0.0005 given for this quasar in Paris et al. (2014). 
We adopt a systematic uncertainty on the redshift of 
±0.0044, or 380 km s“^. This uncertainty is the difference 
between the G lll] emission-line redshift and the principal 
component analysis-based ‘pipeline’ redshift presented in 
Paris et al. (2014); see that reference for details. If our 
adopted redshift is a slight underestimate due to blueshift- 
ing of the emission lines in our spectrum, it is conservative 
in the sense that it errs in the direction of minimizing the 
observed trough outflow velocities. 

J0230 has an apparent magnitude of g — 19.52 and 
an absolute of Mg = —27. Because it is undetected in 
FIRST with an apparent magnitude of i = 18.76, it is not 
radio-loud {Ri < 1; see Figure 19 of Ivezic et al. 2002). 


^ We have determined the features claimed by Foltz 
et al. (1983) to be Ovi at up to 70,000 kms“^ in 
the BAL quasar H 1414-1-089 are actually SlvA1062 and 
S IV*AA1072,1073 absorption in a lower-velocity trough reach¬ 
ing only 28,000 kms“^. That identification is secure because 
the object’s Siv-|-Siv* trough shares the same distinctive 
‘double-dip’ velocity structure as its CIV and N V troughs 
reaching 28,000 kms“^; see Figure 2 of Foltz et al. (1983). 


In section 2 of this paper, we outline our observations, 
data reduction methods, and spectral measurements. In 
section 3 we estimate the mass of the black hole. In section 
4 we compare competing models of BAL variability in the 
context of our mulit-epoch data. Finally, we summarize 
our work in section 5. 

Where needed, we adopt a flat cosmology described 
with Ho = 70 kms“^ Mpc“^, Um=0.3, and 12^=0.7. 

2 SPECTROSCOPIC DATA 
2.1 Observations 

J0230 was identified as having high-velocity absorption 
by visual comparison of its SDSS-I and SDSS-HI spec¬ 
tra as part of a search for newly emerged BAL troughs 
whose results will be reported elsewhere. Follow-up obser¬ 
vations were obtained using the Gemini Observatory (see 
Table 1 for a full list of observations). The Sloan Digital 
Sky Survey (SDSS; York et al. 2000) obtained two spec¬ 
tra of J0230 on MJD 52200 and 52942 using a 2.5 meter 
Ritchey-Ghretien telescope located at the Apache Point 
Observatory in New Mexico. We retrieved the fully re¬ 
duced spectra from the publicly available Data Release 
7 (DR7) quasar catalog (Schneider et al. 2010). These 
spectra cover the wavelength range 3805—9221 A and 
3813—9215 A, respectively, with a spectral resolution of 
R ~ 2000. The Baryon Oscillation Spectroscopic Survey 
(BOSS; Dawson et al. 2013), which is part of SDSS-HI, 
obtained two more spectra of J0230 on MJD 55209 and 
55455 using the same telescope as SDSS, but outfitted 
with a new spectrograph. We retrieved the fully reduced 
spectra from the publicly available Data Release 9 (DR9) 
quasar catalog (Paris et al. 2012). These spectra cover 
the wavelength range 3574—10349 A and 3594—10384 A, 
respectively, with a similar spectral resolution to SDSS. 

We observed J0230 on 15 August 2013 (MJD 56519) 
using the Gemini Multi-Object Spectrograph (CMOS) on 
the 8.1 meter Gemini North Gassegrain telescope located 
at the summit of Manna Kea, Hawai’i. The B600 grating 
with 600 lines mm“^ was used set at a spectrum central 
wavelength of 460 nm (the Blaze wavelength is 461nm). 
Gombined with a 1.0 arcsec wide longslit, the resultant 
wavelength coverage was 3143—6068 A with a spectral 
resolution of i? ~ 1688. The total exposure time was 
1200 s. After noting variability of the broad absorption 
lines in these data, we observed J0230 on 23 December 
2013 (MJD 56649) using CMOS on Gemini North. The 
instrument setup was identical as the previous observa¬ 
tion, but with a longer integration time of 2500 s. The re¬ 
sultant wavelength overage was 3142—6077 A with a spec¬ 
tral resolution of i? ~ 1688. Finally, after noting further 
variability in J0230’s broad absorption lines, we observed 
the object again on 28 January 2014 (MJD 56686) using 
CMOS on Gemini South, the twin telescope to Gemini 
North, located in Gerro Pachon, Ghile. An identical in¬ 
strument setup was used, with a total exposure of 2600 s, 
resulting in a wavelength coverage of 3145—6077 A with a 
spectral resolution of R ~ 1688. All three Gemini obser¬ 
vations were observed at the parallactic angle. The data 
were processed and extracted by standard techniques us¬ 
ing the Gemini IRAF package. The relative fluxes for the 
three Gemini spectra were calibrated using spectropho- 
tometric standard stars; the standard stars were not ob- 
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Observed-frame Wavelength (A) 
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Figure 1. Spectra of J0230 at rest-frame wavelengths (bot¬ 
tom scale) and observed (top scale). The black spectrum is 
the mean SDSS spectra (see § 2.2.1). The blue spectrum was 
taken by BOSS on MJD 55455. The locations of the Civ and 
Si IV emission are labeled (though are weak). In comparing the 
spectra, a broad and deep trough was identified at roughly 
1262—1302 A. This trough was identified as highly blueshifted 
C IV absorption. This trough is referred to as trough A for the 
remainder of the paper. The Flux Density of the BOSS spec¬ 
trum is artificially scaled up to match the continuum level of 
the SDSS spectrum for the purposes of visual comparison. 

served on the same night as J0230. See Table 1 for a list 
of observations and the nomenclature we have adopted 
throughout this paper. 

In Figure 1 we plot the visual comparison of the mean 
SDSS spectrum (black; see § 2.2.1) and the BOSS2 spec¬ 
trum (blue) that led to the identification of an emergent 
absorbing trough. The locations of the C iv and Si iv emis¬ 
sion are given, though noted to be very weak (see § 2.3 
for further discussion regarding Weak Line Quasars). The 
absorbing trough, which we refer to as trough A for the 
remainder of the paper, emerged at some point between 
the two spectral epochs and spans roughly 1260—1300 A. 
We attribute this trough to C iv absorption by highly 
blueshifted gas outflowing along our line of sight to the 
quasar at approximately ~56,000 kms“^. We are con¬ 
fident trough A is not due to blueshifted Siiv absorp¬ 
tion due to the lack of accompanying C iv expected at 
~1425 A. Further, there is some evidence that trough A 
has accompanying N v absorption, which we discuss in 
§ 2.3. We note there is a significant change in J0230’s 
spectrum shortward of trough A, which is attributed to 
changes in the Lya-|-N v complex in that region. Note that 
those changes do not affect our measurements on trough 
A throughout this work. 

2.2 Normalization 

In order to compare the changes we observed in the 
absorption features, we normalized each spectrum by a 
model of the underlying continuum. First, we smoothed 
each spectrum using a boxcar average over 5 pixels. Then 
we identified four regions that appeared unchanged across 
all 7 observations; hereafter, these are referred to as nor¬ 
malization windows. Lastly, we fitted a power-law contin¬ 
uum model to the windows using a least squares routine. 
The normalization windows we used were 1305—1330, 


Observed-frame Wavelength (A) 



Figure 2. The normalized SDSS spectra. While there are some 
small differences between the two spectra (see description in 
§ 2.2.1), since they do not interfere with the two troughs we 
study later we disregard them. As a result, we have combined 
the two SDSS spectra together for the analysis throughout the 
paper. The normalized error spectra are plotted at the bottom. 
The normalization windows are shown as gray regions. 


1410-1420, 1590-1620, and 1650-1700 A. While there 
are other normalization windows that could have been 
used longward of ~ 1700 A for SDSS and BOSS, we re¬ 
frained from using them because the Gemini spectra have 
significantly less spectral coverage. In order to compare 
the SDSS/BOSS spectra to the Gemini spectra we used 
normalization windows accessible from all data. We have 
visually inspected all our spectra and are confident these 
windows represent regions that are unchanged over all 
epochs of our observations. We note this is not a nor¬ 
malization in the traditional sense, as we did not fit the 
emission features. 

Below we describe individual details for normalizing 
the SDSS, BOSS, and Gemini spectra. 


2.2.1 SDSS Normalization 

Two spectra of J0230 were taken on MJD 52200 and 
52942, as part of the SDSS-I survey. In Figure 2 we present 
the normalized SDSS spectra; the gray regions indicated 
the normalization windows. The normalized error spectra 
are plotted along the bottom. Visual comparison show 
little difference between the two, with the small excep¬ 
tion of apparent absorption at 1225—1235 A present in 
the spectrum taken on MJD 52942, but not present in 
the previous epoch, taken on MJD 52200. This feature 
vanished by BOSS2 and never re-appeared, it is present 
only in our noisiest spectrum, and, most importantly, is 
not related to the two broad troughs that are the focus 
of this work. As a result, we do not consider it in this 
study. Other than this feature, there are little differences 
between the two SDSS spectra; we combined them into 
one (hereafter, ’the SDSS spectrum’) in order to increase 
our signal to noise. We adopt an observation date for this 
combined spectrum of MJD 52942, that of the latter SDSS 
observation. Since no broad absorption is present in either 
of the SDSS spectra, we can confidently indicate this date 
to be the last time we observed no absorption present. 
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MJDqjs 

Rest At 

Rest Day 

Plate 

Fiber 

Origin 

SNi676 

Name 

52200.39 

000.00 

- 866.02 

705 

407 

SDSS-I 

7.50 

SDSSl 

52942.34 

213.63 

-652.39 

1509 

365 

SDSS-I 

7.57 

SDSS2 

55208.10 

652.39 

000.00 

3744 

634 

SDSS-III/BOSS 

12.3 

BOSSl 

55454.46 

71.93 

71.93 

4238 

800 

SDSS-III/BOSS 

16.4 

BOSS2 

56519.53 

306.67 

378.06 



Gemini-North 

22.1 

GEMl 

56649.21 

37.33 

415.39 



Gemini-North 

23.2 

GEM2 

56685.07 

10.32 

425.71 



Gemini-South 

18.0 

GEM3 


Table 1. Spectroscopic observations of J0230. Rest At is the rest-frame time in days elapsed since the previous observation. Rest 
Day is cumulative rest days relative to the first BOSS observation. SN 1675 is the median value of the normalized flux divided by the 
error in the flux over the spectral range 1650—1700 A. The final column indicates how we will refer to each epoch for the duration 
of the paper. 


Observed-frame Wavelength (A) 



Figure 3. The normalized BOSS spectra. Since the high- 
velocity absorber at ~1280 A has changed between the two 
epochs, we cannot combine the BOSS spectra. The normalized 
error spectrum for both epochs is plotted at the bottom. The 
normalization windows are shown as gray regions. 


2.2.2 BOSS normalization 

The BOSS survey observed J0230 two more times on 
MJD 55209 and 55455. We normalized these two spec¬ 
tra using the same normalization windows as were used 
for the SDSS spectra. In Figure 3 we plot the normal¬ 
ized BOSS spectra. In both BOSS epochs, trough A is 
present at ~1280 A. The absorption line varies between 
the two BOSS observations, thus we did not combine the 
two spectra as in the case of the SDSS spectra. 


2.2.3 Gemini Normalization 

Three Gemini spectra were taken on MJDs 56519, 56649, 
and 56685. In Figure 4, all three normalized Gemini spec¬ 
tra are plotted. In GEMl we note the emergence of trough 
B, a separate medium-velocity absorber at 1350—1360 A, 
which was not present in any of the SDSS or BOSS spec¬ 
tra. The Gemini spectrum taken on MJD 56685 (orange) 
exhibits less spectral coverage on the red side; the flux 
falls off quickly after ~1675 A. To account for this, the 
third normalization window used for this spectrum was 
1640-1650 A. 


Observed-frame Wavelength (A) 



Figure 4. The normalized Gemini spectra. The normalization 
windows are indicated by the gray regions. The orange spec¬ 
trum (GEMS) has slightly less coverage on the red end and thus 
we changed its third normalization window to 1640—1650 A. 
The normalized error spectra for all three epochs are plotted 
at the bottom. 

2 . 2.4 Final Spectra 

The final six spectra (note the two SDSS spectra were 
combined) are plotted in Figure 5. For reference, the emis¬ 
sion features for Siiv at ~ 1400 A and Civ at ~ 1550 A 
are marked; although both emission lines appear to be 
weak. In our collected data, we note two broad absorp¬ 
tion features, labeled ‘A’ and ‘B’ in the figure. Trough A 
was first observed in the BOSSl spectrum. At its widest 
(BOSS2) trough A spans 40 A (1262—1302 A). Trough 
B was first observed in GEMl. At its widest (GEMl) it 
spans 24 A (1344—1368 A). The legend of Figure 5 in¬ 
dicates the number of rest-frame days since the previous 
observation. 

2.3 Summary of Spectral Features 

In all spectra we obtained of J0230, the emission fea¬ 
tures are relatively weak compared to typical quasars; 
specifically, we measured the rest-frame equivalent widths 
(EWs) of the emission features: Lya-|-Nv= 8.0 ± 0.1 A, 
Siiv < 1.8 A, Civ < 2.5 A, AliiH-Ciii] = 6.1T0.2 A, and 
Mg II = 9.8 ± 0.9 A. When there is no apparent emission 
feature at the expected location of an ion, we measured 
the statistical noise in the spectrum in the ranges provided 
by the Vanden Berk et al. (2001) composite quasar spec- 
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Figure 5. All 6 epochs of spectra plotted together. For ref¬ 
erence, the emission features for Si IV at ~ 1400 A and Civ 
at ~ 1550 A are marked, as well as the two troughs ‘A’ and 
‘B’ we observed to emerge during our monitoring campaign. In 
the legend the MJD of each observation is indicated as well as 
the number of rest-frame days since the previous observation. 
We also note the presence of a third mini-BAL feature labeled 
trough ‘C’ near the systemic redshift of the quasar. There was 
no significant change to trough C through all observations. 

trum (see Table 2 therein). Specifically, they measured 
the Si IV EW over 1360—1446 A, and the Civ EW over 
1494—1620 A. For those regions we measure the statistical 
noise in our spectra to be 0.60 A and 0.84 A, respectively. 
The upper limits quoted above are three times this statis¬ 
tical noise to indicate the largest possible EW these emis¬ 
sion features could have that would still be statistically 
below the noise in our data. Also note that our measure¬ 
ment of Lya-I-N v is contaminated by the Lya forest; the 
actual EW is likely larger. In Luo et al. (2015), a Weak 
Line Quasar (WLQ) is defined as a quasar whose emis¬ 
sion lines have rest-frame equivalent widths of <5 A (they 
drew their sample of WLQs from Plotkin et al. 2010). 
While J0230 does not strictly meet the criterion laid in 
those works, its emission features are still far from a typi¬ 
cal quasar’s. The original WLQ, PG 1407-1-265, has emis¬ 
sion features with comparable EWs to J0230 (McDowell 
et al. 1995, see Table 2 therein), as does the prototypical 
WLQ, PHL 1811, which which has the following EWs: 
Lya-I-Nv= 15 A, Civ = 6.6 A, AliiH-Ciii] <4 A, and 
MgII = 12.9 A (Leighly et al. 2007). Further, quasars 
with EWs < 10 A investigated so far have sufficient sim¬ 
ilarities (e.g., common X-ray weakness) and can likely be 
unified as per Luo et al. (2015) in a common physical 
model. Therefore, we consider J0230 a WLQ.® 

Also present in all spectra is a narrow C iv absorption 
feature at ~1550 A, very close to the systemic redshift of 
J0230 (also seen in Siiv, Cll, Nv, and Lya), hereafter 
trough C. There were no dramatic changes in trough C in 
our observations. 

® WLQs tend to have blueshifted broad emission lines in the 
UV, making systemic redshift determination more challeng¬ 
ing than usual. Our adopted systematic redshift uncertainty of 
±380 kms“^ in J0230 is similar to the ±300 kms“^average 
difference between redshifts determined by narrow-line studies 
and those determined by SDSS for weak line quasars found by 
Plotkin et al. (2010). 


The changes in the spectrum are best seen in Fig¬ 
ure 6. In BOSSl we note the appearance of trough A: 
a broad, high-velocity absorber covering the wavelength 
range 1260—1300 A. Trough A grew to its strongest in 
BOSS2 by getting both deeper and wider; these changes 
were mostly in the the low-velocity half of the trough, 
whereas the high-velocity half of the trough changed 
less. In the first Gemini spectrum (GEMl), the high- 
velocity half of trough A weakened greatly while its low- 
velocity half weakened only somewhat, in comparison to 
BOSS2. Between GEMl and GEM2, trough A strength¬ 
ened slightly on its high-velocity side. We note in GEMl 
the emergence of trough B, a second high-velocity ab¬ 
sorber in the wavelength range 1344—1368 A. We are 
confident this absorption is due to highly blueshifted 
gas along the line of sight to J0230. It cannot be due 
to blueshifted Silv absorption because that would re¬ 
quire accompanying Giv at ~1500 A. Further, there is 
some evidence to suggest there is accompanying Si iv and 
N v at similar outflow velocities (see below)y. Trough B’s 
low-velocity end remained relatively unchanged (though 
slightly weaker) into GEM2, while its high-velocity side 
reached higher outflow velocities. Finally, between GEM2 
and GEM3, trough A did not change appreciably, while 
trough B weakened on its low-velocity side and its high- 
velocity edge reached higher outflow velocities. Trough B 
also decreased in depth. 

The presence of C iv absorption can be accompanied 
by absorption of one or more other ionic transitions, such 
as Si IV, Lya, and N v. We searched for absorption of these 
ions that would correspond to the same outflow velocities 
as trough A or B. Figure 7 shows all 6 normalized spectra 
with a much heavier smoothing, and with a much wider 
wavelength coverage. We have marked the observed loca¬ 
tions of the G IV absorption by trough A (dashed line) and 
by trough B (sold line), along with the expected locations 
of their accompanying Silv, Lya, and Nv absorption. We 
have plotted the error spectra of the SDSS, BOSS, and 
GEMl spectra along the bottom. For the purposes of clar¬ 
ity, the spectra were heavily smoothed in order to see fea¬ 
tures better in this more noisy part of the spectra. It is 
also of note the normalization was not repeated with new 
normalization windows in the region from 1000—1300 A, 
thus the relative flux levels are not necessarily accurate. 
This is only meant to be a search for possible accompa¬ 
nying absorption. 

For trough A, there appears to be no accompanying 
Si IV absorption in any of the spectral epochs we obtained. 
In searching for accompanying Lya±N v, we note that the 
wavelength coverage does not extend far enough into the 
blue for SDSS, BOSSl, or BOSS2 but does for the three 
Gemini spectra. In these latter three epochs there may be 
N V, but no apparent Lya is observed. For trough B, we 
note the possible presence of accompanying Silv absorp¬ 
tion in the three Gemini spectra, however, the absorp¬ 
tion is coincident with the Lya±N v emission systemic 
to J0230. Since it is very difficult to disentangle emission 
from coincident absorption, we cannot confirm this to be 
Si IV. The identification is also not certain because a flux 
deficit was also seen at that location in the SDSS spec¬ 
trum, before trough B appeared. There is probable N v 
absorption for trough B. 

Archival photometry of J0230 is available since it is 
located in Stripe 82, a region of sky imaged by SDSS, mul- 
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Figure 6. Each spectral epoch is plotted centred on the two 
absorbers, troughs A and B. We have separated the spectra ar¬ 
tificially by 0.5 normalized flux units, with the earliest epoch 
(SDSS) at the top, and the most recent (GEMS) at the bot¬ 
tom. The dashed lines indicate the normalized continuum level 
for each spectrum. The black bars indicate where we define the 
absorption features to begin and end. Note for the SDSS spec¬ 
trum, there is no apparent absorption in either troughs A or 
B. Also note for SDSS, BOSSl, and BOSS2 there is no appar¬ 
ent absorption for trough B. For these cases we have placed 
a slightly thinner black bar across the regions that represent 
the widest that trough became. For trough A, this occurs in 
BOSS2 and for trough B this occurs in GEMl. 

tiple times over 7 years (MacLeod et al. 2010). We have 
obtained the photometry of J0230 from the SDSS archive, 
however, it is not concurrent with our spectroscopy. Thus 
it cannot help ns interpret the spectroscopic variability we 
observe. J0230 was too faint for the Catalina Real-time 
Transient Survey (CRTS). 

2.4 Measurements of Troughs A and B 

We measure the properties of absorption troughs A, B, 
and C, such as the equivalent width (EW), the weighted 
centroid velocity Ucent, and the average trough depth, in 
all observations in order to compare changes from one 
epoch to the next. 

In Figure 6, the 6 epochs of normalized spectra are 


plotted (separated artificially in the y direction). The bot¬ 
tom x-axis is the rest-frame wavelength, and the top x- 
axis is the outflow velocity relative to C iv ~1550 A. The 
dashed lines indicate the continuum for each spectrum. 
The dark horizontal lines indicate where we define ab¬ 
sorption is present for troughs A and B (see below on 
how these were chosen). 

To measure the EW from the normalized spectra we 
followed equations 1 and 2 in Kaspi et al. (2002), which 
are, 

EW = E(l-f)^- (1) 

and the uncertainty on the EW is. 



Fi and AT) are the flux and the error on the flux in the 
ith bin, is the underlying continuum flux, AT), is the 
uncertainty in the mean of the continuum flux in the nor¬ 
malization windows, and Bi is the bin width in units of 
A. In our normalized spectra, F^ = 1 and AEc is calcu¬ 
lated using the windows 1305—1330 A and 1410—1420 A, 
which are the two windows closest to both absorption 
troughs. Thus ctew represents the statistical uncertainty 
inherent to the spectra. It does not quantify the systemic 
uncertainty, which is governed by the placement of the 
continuum by normalization. 

The edges of troughs A and B in a given spectrum 
were identified by finding the locations where the flux 
drops below, and stays below, the normalized continuum 
level of Ec = 1. In Figure 6, these edges are represented by 
black horizontal bars; in Table 2, AW is calculated using 
these edges. We applied Equations 1 and 2 to calculate 
the EW within the edges found. We note that the place¬ 
ment of the normalized continuum, and thus the locations 
of the edges of the troughs, is highly sensitive to the nor¬ 
malization process. Further, for trough A, the absorption 
appears to be truncated by the Lya-I-N v emission com¬ 
plex; as a result we consider our EW measurements to be 
conservative. 

Note for both troughs, some epochs do not exhibit 
absorption; both troughs A and B are not present in the 
SDSS spectrum, and trough B does not exhibit absorp¬ 
tion in the SDSS, BOSSl, and BOSS2 spectra. For these 
cases, we took the largest trough width determined for 
each trough and applied it to the unabsorbed spectra. 
For example, in the case of trough A, the widest the 
trough was observed to be was in the BOSS2 spectrum 
at 1262 — 1302 = 40 A. We applied this range of the ab¬ 
sorption profile in the unabsorbed spectra of SDSSl and 
measured the EW. The resulting value for the SDSS spec¬ 
trum was —0.18±0.48, indicating an EW consistent with 
zero. More examples of this can be found in Table 2 la¬ 
beled with an ellipsis. 

We measured the centroid velocity, Vcent, of the 
trough following the definition in Filiz Ak et al. (2013); it 
is the mean of the velocity in a trough where each pixel is 
weighted by its distance from the normalized continuum. 

The mean depth of the trough was calculated in two 
ways. First, we measured (Ibal as in Filiz Ak et al. (2013), 
which is the mean distance from the normalized contin¬ 
uum level for each data point in the trough. Second, we 
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Figure 7. Plotted are the 6 normalized spectra with heavier smoothing (boxcar with window of 25 pixels) and with a much 
wider wavelength coverage. We have marked the location of the trough A C IV absorption with a vertical dashed line; the expected 
locations of Si IV, N V, and Lyoabsorption features that may accompany trough A’s C IV are also marked with vertical dashed lines. 
The location of trough B’s C IV absorption is marked with a vertical solid line, as are the expected locations of this features possible 
accompanying Si IV, N V, and hya. At the shortest wavelengths of the BOSSl and BOSS2 spectra (at <1330 A in the trough A 
rest-frame), the spurious broad ’emission’ features are due to noise. It is of note these are the normalized spectra from § 2.2, which 
were created using continuum windows between 1300—1700 A. 


measured dmaxv, which is calculated by sliding a 7 pixel¬ 
wide window across the trough and measuring the aver¬ 
age depth over each window. We take the largest value 
of all these windows as d,nax 7 - The uncertainty on the 
depth calculated as the uncertainty in the mean of the 7 
pixels in the average. We note that since the our observa¬ 
tions were taken with different telescope and instrument 
set ups, 7 pixels corresponds to slightly different resolu¬ 
tions; however, the differences are too small to impact the 
measurements. For reference, the 7 pixels covers approx¬ 
imately 2 A, or ~450 kms“^ in all spectra. 

2.4-1 Coordinated Variability 

Work on BAL quasar variability indicates troughs from 
the same object can vary in coordination with each other, 
which can lead to constraints on variability models (i.e., 
Filiz Ak et al. 2012, Wang et al. 2015; see discussion in 
§4.1.1 below). We plot the EW of each trough vs. the rest- 
frame time elapsed since the SDSS epoch in Figure 8 and 
dmax 7 for each trough vs. the rest-frame time elapsed since 
the SDSS epoch in Figure 9 in order to investigate how 
the variability of one trough compares with the others. 
The EW of trough C remains relatively constant over all 
epochs. Both trough A and B begin with a very low EW, 
then emerge with a sharp and significant increase in later 
epochs (BOSSl for A and GEMl for B). There is an inter¬ 
esting pattern in the final three observations (the Gemini 
epochs), which occurs after both troughs have emerged 
and are established: the EW for both trough A and B 


increases for GEM2 and then returns to the same value 
it was in GEMl for GEMS. This pattern could be inter¬ 
preted as absorption from two physically distinct clouds 
varying in a coordinated fashion (for reference, the time 
frame from GEMl to GEMS is 47 days). However, the un¬ 
certainties on our EWs are of similar scale to the amount 
of variability we are referring to in the Gemini epochs. 
Thus, this pattern does not represent statistically signifi¬ 
cant coordination in variability. 

2-4.2 BALnicity Index 

For comparison to other BAL-quasar studies, we have 
measured the BALnicity Index (BI) of J02S0. We calcu¬ 
late the Absorption Index (AHso), defined in Hall et al. 
(2002), following: 



where f{v) is the normalized flux density as a function 
of velocity, and C' is equal to 1.0 within a trough if 
the trough is wider wider widerthan 450 kms~^, oth¬ 
erwise it is set to 0.0. The integration begins at n = 
0 kms~^ relative to the systemic velocity of the quasar 
and runs through the highest velocity at which absorption 
is present. 

We also measure the modified BALnicity index, BI*, 
following: 
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Figure 8 . The measured equivalent widths (EW) for troughs 
A (blue), B (green), and C (red). The trough C points are 
artificially shifted to the right by 10 days in order to avoid 
confusion with trough A data points. 


Figure 10. Combined BOSS spectra (black), continuum fit 
(dashed black), and Gaussian fit (dashed red) to the Mg II emis¬ 
sion feature at ~ 2800 A. The fit was applied only to the data 
in grayed out region. The best-fit Gaussian parameters to the 
data are shown in the lower left. 
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Figure 9. The maximum depth of trough A (blue), B (green), 
and G (red) as a function of rest-frame days. dmax 7 represents 
the lowest average 7 pixels in a row for each trough. 


which is a modification of the original BI defined in Wey- 
mann et al. (1991) that imposes no constraint on the out¬ 
flow velocity. BI* has no formal limit on the minimum and 
maximum absorbing velocities. The quantity C is equal 
to 1 only when the quantity in parentheses is greater than 
zero for more than 2000 kms“^, otherwise it is set to 0.0. 

In Table 3, we list the BALnicity indexes calculated 
using both methods. The total index value is measured 
over viow > 0 and Vhigh < 60,500 km s“^, however, we 
also provide the individual contributions of each trough 
in the spectra. Note that for AI450, trough C contributes 
to the total index, but for BI* it does not. 


3 BLACK HOLE MASS ESTIMATE 

To estimate the mass of the SMBH in J0230, we used the 
velocity dispersion of the Mg 11 A 2796, 2803 A emission 
line. A full description of this technique can be found in 
Rafiee & Hall (2011). Equation 9 of that work is 

Mbh/Mq = 30.5[AL30oo/( 10^^ erg s-^)]'’ ®cr^ (5) 


where Laooo is the observed monochromatic luminosity at 
3000 A rest-frame, A = 3000 A, and a is the intrinsic 
line dispersion of the Mgll emission line in kms~^. There 
is intrinsic scatter of ±0.15 dex (±35%) and systematic 
uncertainty of ±0.10 dex (±24%) in this equation. 

The two BOSS spectra of J0230 represent the best 
coverage we have of that wavelength regime. We combined 
the two BOSS spectra with a weighted mean, and then 
fit a line to the continuum using windows 2650—2700 A, 
2900—3000 A. After fitting and removing the continuum, 
we fit a Gaussian to the remaining Mgll emission in the 
region 2700—2900 A. In Figure 10 the fitted Gaussian 
is plotted over the normalized BOSS spectra. The best- 
fit parameters were fi = 2805 A, and a = 22.1 A. The 
standard deviation in the Gaussian fit indicates the ve¬ 
locity dispersion of the Mgll emission feature, which is 
caused by the Doppler broadening of an AGN broad line 
region orbiting the SMBH. We convert cr = 22.1 A = 
2370 kms"A 

To calculate the quasar luminosity we used 

^L^ooo = dTrD^/aooo x 3000(1 ± z), (6) 

where Dl is the luminosity distance, /3000 is the observed 
flux density at rest-frame 3000 A, and z is the redshift. We 
measured /3000 from the combined BOSS spectrum to be 
/3000 = 3.35x 10“^^ erg s“^ cm“^ (observed A)“^. For our 
adopted cosmology, the luminosity distance to z = 2.473 
is Dl = 2.01 X 10^° pc, or 6.21 x 10^® cm. Therefore, we 
have AL3000 = 1-69 x 10^® erg s“^. Finally, we calculate 
the mass of the SMBH to be Mbh/Mq = 2.2 x 10® with 
an intrinsic scatter of ±0.8 x 10® (±35%). 

For a SMBH of this mass the Eddington Luminos¬ 
ity is LEdd = 3.45 X 10“^^ erg s“^. Using a bolomet- 
ric correction of BC3000 = 5 derived in Richards et al. 
(2006), this quasar has an estimated LboI = 8.45 x 10"^® 
erg s“^, and therefore this quasar has an estimated 
fsdd = LBoi/LEdd = 0.25. Such a black hole has Rsch = 
6.6 X 10® km. 

It is worth pointing out that Plotkin et al. (2010) 
present some evidence to suggest that some Mgll emis¬ 
sion lines of WLQs could exhibit non-virialized behaviour 
(namely, the emission feature is blueshifted from the sys- 
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Rest At 

EWibcr 

Aw 

Av 

dmax7 

Hcent 

dsAL 

Trough A 

(days) 

(A) 

(A) 

(km/s) 


(km/s) 


SDSS 

000.00 

-0.19±0.48 



0.07±0.06 



BOSSl 

652.S9 

3.23±0.34 

31 

7002 

0.18±0.05 

56496 

0.10 

BOSS2 

71.9S 

6.68±0.30 

40 

9028 

0.26±0.04 

56004 

0.17 

GEMl 

306.67 

2.72±0.41 

27 

6063 

0.18±0.08 

53769 

0.10 

GEM2 

37.33 

3.45±0.31 

26 

5860 

0.21±0.06 

55020 

0.13 

GEM3 

10.32 

2.70±0.28 

23 

5185 

0.22±0.06 

55101 

0.12 

Trough B 

SDSS 

000.00 

0.92±0.31 



0.06±0.05 



BOSSl 

652.39 

0.37±0.28 



0.03±0.05 



BOSS2 

71.93 

0 . 21 ± 0.22 



0.03±0.04 



GEMl 

306.67 

3.33±0.31 

24 

5214 

0.21±0.05 

39212 

0.14 

GEM2 

37.33 

4.26±0.25 

22 

4780 

0.31±0.05 

39726 

0.11 

GEMS 

10.32 

3.59±0.21 

20 

4352 

0.27±0.04 

40224 

0.10 

Trough G 

SDSS 

000.00 

2.89±0.13 

8 

1549 

0.49±0.04 

87 

0.37±0.05 

BOSSl 

652.39 

2.67±0.12 

8 

1550 

0.47±0.04 

78 

0.33±0.05 

BOSS2 

71.93 

2.98±0.10 

10 

1935 

0.52±0.03 

68 

0.30±0.05 

GEMl 

306.67 

2.74±0.08 

7 

1356 

0.53±0.03 

163 

0.39±0.04 

GEM2 

37.33 

2.97±0.07 

7 

1356 

0.56±0.02 

125 

0.42±0.04 

GEMS 

10.32 

2.74±0.08 

7 

1356 

0.53±0.03 

212 

0.39±0.04 


Table 2. Measurements made on trough A, B, and C. The ’ indicate where no absorption is visible in the spectrum. Values of 
EW for these cases used the widest possible AW the trough was observed to reach (BOSS2 for trough A, GEMl for trough B). 



AIa 

AIs 

Ale 

total Al 

BIl 

BI^ 


total BI* 

SDSS 

0.0 

0.0 

477±3 

477±3 

0.0 

0.0 

0.0 

0.0 

BOSSl 

152±5 

0.0 

437±2 

589±5 

0.0 

0.0 

0.0 

0.0 

BOSS2 

746±6 

0.0 

490±2 

1236±6 

561±5 

0.0 

0.0 

561±5 

GEMl 

103±6 

323±6 

455±2 

880±8 

0.0 

115±4 

0.0 

115±4 

GEM2 

293±6 

547±6 

491 ±2 

1331±9 

74±4 

242±4 

0.0 

316±5 

GEMS 

210±6 

433±5 

453±2 

1096±7 

35±3 

155±3 

0.0 

190±4 


Table 3. The BALnicity was calculated using two different definitions: AI 450 and BI* (see § 2.4.2). We calculated the total index 
over a velocity range of > 0 and < 60,500 km s“^. We also calculated the individual contributions to the index by 

each trough in the spectra. In the case of AI 450 , trough C contributed to the measurement; for completeness, we provide its index 
measurement. In BI*, trough C did not contribute to the total. Note the uncertainties quoted here are statistical only. Systematic 
uncertainty introduced by the placement of the continuum is not taken into account. A reasonable continuum uncertainty of ±5% 
translates to a BALnicity Index uncertainty of ±5%ldBAL- 


temic redshift, though only by 360 kms”'^ on average; see 
§ 5.2 and 6.1 of that work). We see no such evidence of 
a non-virialized Mgll emission feature in J0230: a single 
Gaussian function fits the emission line well, its peak is 
actually redshifted by ~ 510 ± 380 kms“^ from the po¬ 
sition of Mg II expected from the composite spectrum of 
Vanden Berk et al. (2001). 


Moreover, we have calculated black hole masses us¬ 
ing the dispersion in the Mgll emission line for the ob¬ 
jects from Plotkin et al. (2010). We find that the resulting 
masses are larger than the black hole masses they calcu¬ 
late using the dispersion in H /3 by only a factor of two. 
A deviation of that factor is not statistically significant 
given the uncertainties on our black hole mass estimate. 


4 DISCUSSION 

As mentioned in the introduction, broad absorption 
trough variability in quasars can be explained by trans¬ 
verse motion of the absorbing clouds across the line of 
sight to the accretion disk, or by changes in the ionization 
parameter of the absorbing cloud, or by a combination of 
these. Here we analyze two possibilities individually lay¬ 
ing out constraints where possible. Note that in the end, 
the range of possible locations for the gas is large enough 
to preclude useful constraints on the kinetic luminosity of 
the outflow (Dunn et al. 2010), especially since the solid 
angle covered by this extremely high velocity outflow is 
unknown. 


4.1 Pure Transverse Motion Variability Model 

In the transverse-motion model it is assumed the absorp¬ 
tion parameters of the cloud of gas are unchanged, and all 
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changes to the EW, the velocity profile, and the maximum 
depth of the trough can be explained by an absorbing 
cloud moving to cover more or less of the accretion disk. 
Any evolution of an absorption feature (i.e., an emergent 
trough) can be explained as long as the constraints from 
timescales yield plausible transverse velocities. 

The transverse velocity of an absorbing cloud across 
the line of sight is derived by dividing the distance the 
cloud travels by the travel time it took to get there, the 
latter of which is time between successive observations. 
In order to measure the distance covered by an absorbing 
cloud between those observations we must both estimate 
the size of the continuum region it is traversing, and also 
model the relative sizes and shapes of the cloud and con¬ 
tinuum region. 

We approximate the continuum region to be repre¬ 
sented by the a-disk model presented in Shakura & Sun- 
yaev (1973), hereafter SS73, with the following model pa¬ 
rameters. We set a = 0.1, a free parameter in the model 
that governs the amount of accretion as a result of tur¬ 
bulence, typically 0 < a < 1. We assume a non-spinning 
black hole, which leads to an accretion efficiency of rj = 
0.0572. Given these parameters, the rate of mass accre¬ 
tion onto the black hole would be rh = fEdd/v = 0.25. 
Using an accretion disk defined by these parameters, we 
can estimate 1195(1320), the continuum diameter within 
which 95% of the 1320 A continuum is emitted. We use 
the 1320 A continuum, which is the region in between 
troughs A and B, because it allows us to use a conin- 
tuum region that is the same size for both troughs; we 
note the size of the accretion disk would only change 
a small percentage if using the trough A or B centroid 
wavelengths. We find 1)95(1320) = 63 Rsch, therefore, 
7 ) 95 ( 1320 ) = 4.2 X 10^^ km. That gives a light-crossing 
time of 1.4 X 10® s = 16 days. 

However, accretion-disk sizes inferred from 
gravitational-microlensing studies and photometric- 
reverberation studies (e.g., Morgan et al. 2010, Black- 
burne et al. 2011, Jimenez-Vicente et al. 2012, Edelson 
et al. 2015) are approximately a factor of four larger than 
the theoretical size predicted in the SS73 a-disk model 
(see a full discussion in Hall et al. 2013). Therefore, we 
increase our estimated continuum-source diameter by a 
factor of four, to 7 ) 95 ( 1320 ) = 252 Rsch = 1-7 x 10^^ km. 
The uncertainty in this number is likely a factor of two. 
A disk that size has a light-crossing time of 5 x 10® s = 
64 days. 

With the estimated size of the emitting region, and, 
given some simple models of clouds moving into or out of 
the line of sight of an emitting region, we can estimate a 
maximum and minimum transverse velocity of an absorp¬ 
tion cloud that would be responsible for the emergence 
and variability of troughs A and B. 

The most dramatic change we observed in the absorp¬ 
tion depth of J0230 occurred in trough B when it emerged 
between the BOSS2 and GEMl observations; the change 
in depth was Admax? = 0.21 — 0.03 = 0.18 over a period of 
307 rest-frame days. As per the transverse-motion model, 
if we consider this change in depth to be entirely due to 
more of an optically thick absorbing cloud moving into 
the line of sight to the emitting region, it suggests that 
over 307 rest-frame days, the emitting region went from 
having 3% of its flux blocked to having 21% blocked, or 
a 21% covering fraction, C. Note in order for changes in 


absorption depth to equate to changes in covering frac¬ 
tion we are assuming the lines are optically thick. (If the 
lines are optically thin, the absorber must reach a larger 
covering fraction of the emission region in the same time 
span, requiring even higher transverse velocities.) 

In Gapellupo et al. (2013), two simple models were 
proposed for clouds crossing the emitting region (see Fig¬ 
ure 14 therein). The first scenario is the ‘crossing disks’ 
model, where the absorbing cloud is projected on the sky 
as a circle (or a disk) and is crossing a circular emitting 
region (where the emitting region appears much larger 
than the absorbing cloud). In the second scenario the 
absorber is much larger than the background emitter 
it is traversing; this is the ‘knife-edge’ model. As men¬ 
tioned above, the crossing speeds in these two scenar¬ 
ios are measured by dividing the distance traveled by 
the time it took to travel there. The change in cover¬ 
ing fraction, AC, is the fraction of the emitting region 
the absorber crosses in the time-frame between obser¬ 
vations. Therefore, in the ‘crossing disks’ scenario, the 
minimum distance traveled by the gas responsible for the 
emergence of trough B is x/AG 7)95(1320), and the cross¬ 
ing time is At = 307 days. Therefore, VO.18 x 64 light- 
days = 27 light-days in 307 days. Therefore the trans¬ 
verse speed is 26,400 kms“^. However, if we assume the 
cloud has traversed to the exact opposite side of the emit¬ 
ting region, the distance traveled is the complete 64 light- 
days in 307 days. This results in a transverse velocity of 
62,500 kms“^. In the ‘knife-edge’ scenario, the distance 
traveled is AG 7 ) 95 ( 1320 ) = 12 light-days in 307 days. 
This equates to 11,700 kms”*^. Thus, given the above two 
scenarios, we can place the transverse velocity of a cloud 
responsible for the emergence of trough B in the range 
11,700 < u(kms“^) < 62,500. For trough A, the most 
dramatic change in absorption depth also occurred dur¬ 
ing its emergence, which was between SDSS and BOSSl; 
the change in depth was AdmaxT ~ 0.18 — 0.07 = 0.11 
over a period of 652 rest-frame days. Applying the same 
relations as above we can place the transverse velocity of 
a cloud responsible for the emergence of trough A in the 
range 3, 200 < u(kms“^) < 29, 500. 

While these two models can be useful in interpreting 
observations in a campaign with two epochs, our unique 
dataset consists of six epochs. Analyzing the behaviour 
of the absorption features over all six epochs allows us to 
test the predictive power of the above two scenarios. For 
instance, trough B was consistent with zero absorption in 
the SDSS, BOSSl, and BOSS2 observations (see Table 2). 
The trough appeared between the observation of BOSS2 
and GEMl, which was over a time period of 307 days, 
then for the next 2 observations (GEM2, and GEM3) the 
trough remained close to the same depth and EW (within 
the uncertainties). Assuming an absorber is moving at a 
constant velocity transverse to the line of sight, the above 
behaviour rules out the ‘knife-edge’ scenario, which would 
only cover more area as time goes on. 

If we assume the emitting region has a uniform flux 
across its area (as Gapellupo et al. 2013 does), then the 
‘crossing disk’ scenario can explain trough B’s behaviour. 
However, research into the theoretical understanding of 
accretion disks - through the work of SS73 and Dexter & 
Agol (2011) (among others) - indicates the emitting region 
is unlikely to be homogeneous. If we assume the emitting 
region is more like a SS73 disk, where the majority of the 
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flux is concentrated toward the centre of the emitting re¬ 
gion, we can also rule out the ‘crossing disks’ scenario.^ 
A crossing disk of fixed size traversing a SS73 accretion 
disk at a constant velocity would produce an increasing 
amount of observed covering fraction as it moved across 
the first half of the disk, but then a decreasing amount 
of covering as it traversed the second half of the disk. 
If trough B appeared in GEMl as a result of transverse 
motion, we would have expected to see the depth of the 
absorber decrease appreciably in the subsequent obser¬ 
vations of GEM2 and GEMS. Since this is not the case, 
the ‘crossing disks’ scenario is unlikely to be the correct 
interpretation of the variability of trough B. 

Over 6 epochs, the nature of trough A’s variability 
also rules out the ‘knife-edge,’ but agrees with the aug¬ 
mented ‘crossing disks’-|-SS73 scenario. Specifically, there 
was no measured absorption in SDSS after which there 
was an increase in absorption in BOSSl which continued 
to increase in both depth and EW into BOSS2. Then by 
GEMl through GEMS, both the depth and EW returned 
back to values similar to those measured in BOSSl. This 
is consistent with a cloud smaller in angular size than 
the emitting source traversing into the line of sight for 
BOSSl, crossing the central portion of the disk leading 
to a the measurements of BOSS2, continuing on to the 
second half of the disk for GEMl and GEM2, but has 
not reached the other side yet as there is still measured 
absorption in GEM3. If we apply the relations from the 
‘crossing disks’ scenario above, the velocity range this 
absorber would have is 10,000 < v (kms“^) < 18,000. 
At 18,000 kms“^ we expect trough A to disappear com¬ 
pletely approximately 350 days after the GEM3 observa¬ 
tion. At 10,000 kms“^ we expect it to disappear approx¬ 
imately 1,500 days after the GEM3 observation. 

There is one plausible scenario of transverse motion 
that can match the observed depth changes in trough B 
in the context of a SS73 disk: a ‘flow-tube’ (similar to 
that proposed by Arav et al. 1999; see their Figure 10). 
In Figure 11, we have plotted a log-luminosity map of an 
accretion-disk emitting at 1320 A powered by a SMBH 
equal to that of J0230 (see § 3). The emitted light is much 
more concentrated towards the centre (though note there 
is a region occupied by the black hole where no emis¬ 
sion is observed). We have plotted over top of the map 
an example of our proposed flow-tube scenario. The tube 
is traversing the continuum region at some impact pa¬ 
rameter, i, away from the centre, and has some width, 
w. The tube extends infinitely to the left in this figure. 
We note that our flow-tube geometry and dynamics differ 
from that proposed in Arav et al. (1999). Specifically, we 
have chosen a flow-tube that is homogeneous from centre 
to edge and is in the midst of establishing itself along our 
sight-line before settling into a long-term configuration as 
discussed in Arav et al. (1999). 

If a flow-tube similar to the one shown in Figure 11 
were to move across the emitting region of J0230, it would 
serve to create a sharp increase in absorption as it crossed 
close to the centre of the disk, but due to there being very 
little flux at the edges of a SS73 disk, not much more 
coverage would occur as it traversed the second half of 

^ see Fig. 11 for an example of the luminosity gradient of a 
SS73 accretion-disk; this figure will be discussed in more detail 
later 
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Figure 11. An example of a flow-tube traversing a simulated 
emitting region of an accretion-disk. The logarithm of the lumi¬ 
nosity of the disk is represented by the gray scale. Over-plotted 
is an example of a flow-tube traversing the disk, which would 
serve to cover some of the light, creating absorption. In this 
representation, the tube extends infinitely to the left but ter¬ 
minates at the right edge drawn. The width of the tube is w 
and the impact parameter relative to the centre of the accre¬ 
tion disk is i. Note there is a region at the very centre occupied 
by the black hole where no luminosity is observed. 


the continuum region. This geometry would match the 
behaviour we see in the variability of trough B. 

We have investigated whether a flow-tube of this na¬ 
ture could successfully reproduce the variability in trough 
B, and at what velocities it could do this, by simulat¬ 
ing flow-tubes of various widths and impact parameters 
traversing a SS73 disk, measuring how much flux is cov¬ 
ered as a function of distance across the disk the simu¬ 
lated flow-tubes produce, and then attempting to match 
the observed covering fractions for trough B to the simu¬ 
lated covering fraction vs. distance generated by the flow- 
tubes. Referring to Figure 11, traversed distance is mea¬ 
sured along the x-axis of the disk, and the direction of 
motion of each simulated flow-tube is from the negative 
X direction towards the positive x direction. 

Matching the observations to our simulations was 
done via the following prescription: A given simulated 
flow-tube has covering fraction as a function of x, C{x). 
We search for a distance across the accretion disk, xo, 
that matches the covering fraction for BOSS2, namely 0 < 
C{xo) < C(BOSS2)-|-l(j, which is the last time trough B 
was measured to have a depth consistent with zero. When 
found, we go searching for the next closest x\ that sat¬ 
isfies C{GEM1) - lo- < C( 2 ;l) < C{GEM1) + la. We 
calculated the velocity, v, required to cover the distance 
from Xo to xi, given the known time between successive 
observation (307 days from BOSS2 to GEMl). We then 
searched for the next X 2 that satisfies C{GEM2) — Ict < 
G{x 2 ) < G{GEM2) -b la. When a match is found, we use 
the simulated distance from xi (GEMl) to X 2 (GEM2), 
and the velocity the flow-tube is moving at, v calcu¬ 
lated above, to determine the length of time it would 
take for the flow-tube to cover the distance xi — X 2 . If 
the time is equal to the time between GEMl and GEM2 
observations (37 days) then we continue the search to 
see if GEM3 also matches. We look for 0:3 that satis¬ 
fies G{GEM3) - la < G{x 3 ) < G{GEM3) -b la. Sim- 
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impact parameter (R.s,,,) 


Figure 12. Width of flow-tube vs. the distance from centre of 
accretion disk the flow-tube traverses (impact parameter). It is 
plotted in units of Rsch = 6.6 X 10® km. The gray region repre¬ 
sents all possible combinations of flow-tubes that resulted in a 
final covering fraction between 15% and 30%. The black points 
are the combinations of parameters that not only matched all 
covering fractions in our observations, but did so within the 
observation time constraints. The x-axis is plotted is distance 
from centre of tube to centre of accretion disk, where positive 
and negative values represent opposite sides of centre. 

ilar to above, we use the distance from X 2 (GEM2) to xs 
(GEMS), and the v above to determine the length of time 
it would take for the flow-tube to cover that distance. If 
that time is equal to the time between GEM2 and GEMS 
observations (10 days) then we have found a combina¬ 
tion of width and impact parameter for a simulated flow- 
tube that matches the variability in covering fraction as 
well as the time between successive observations. In Fig¬ 
ure 12, we have plotted the parameter space of width vs. 
impact parameter that we investigated with the simulated 
flow-tubes. The gray region displays the combinations of 
parameters that resulted in a flow-tube’s final covering 
fraction (after it had completely traversed the disk) be¬ 
tween 15% and S0%, which is a healthy margin for the 
GEMS covering fraction. The black points represent the 
combinations that fit the variability of BOSS2 through 
to GEMS. In Figure 13 we plot a histogram of all possi¬ 
ble velocities we determine from the above analysis. The 
mean velocity of the distribution is 36, 800 kms“^ with a 
range spanning 8, 000 < v (kms“^) < 56, 000. 

Trough A could still be explained as a flow tube, but 
a simple flow tube model is not consistent with its dmax 7 
and EW variations. The best fit despite those variations 
would yield a slower transverse velocity because the time 
over which the biggest change occurred (SDSS-BOSSl) is 
larger than for trough B. Note that a slower transverse 
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Figure 13. The range of possible velocities of a flow-tube 
traversing the emitting region of J0230. These were determined 
by simulating flow tubes of various widths and impact param¬ 
eters across a SS73 disk scaled to match J0230’s mass and 
monochromatic luminosity at 1320 A. In order to be a plau¬ 
sible velocity, the tube must recreate the covering fraction at 
each spectral epoch, given one velocity, as well as match the 
time between observations. 

velocity is consistent with trough A’s higher line-of-sight 
velocity, as gas which is closer to terminal velocity is likely 
farther from the quasar with lower transverse velocity due 
to angular momentum conservation. 

In summary, we have found that pure transverse mo¬ 
tion can plausibly explain the variability of both trough 
A and B over all 6 epochs of observation in our dataset. 
Trough A is best explained by a ‘crossing disks’ traversing 
a SS73 disk at velocities between 10,000 < v (kms“^) < 
18, 000. This model and velocity range allow us to predict 
trough A will disappear between 350 < t (days) < 1,500 
after our last observation (GEM3). Trough B is best ex¬ 
plained by a flow-tube that has recently moved into the 
line of sight, travelling in the velocity range 8,000 < 
V (kms”*^) < 56,000. In this scenario, we have no con¬ 
straint on how far a flow-tube extends, and thus cannot 
predict when trough B will disappear. 



4 . 1.1 Constraining Distances 

If we assume the absorbers responsible for both trough A 
and B have reached maximum velocity, have transverse 
velocities vtrans and we are observing them at some cur¬ 
rent distance rc from the BH, but was launched from 
a circular orbit at a distance ri,, which has an orbital 
velocity {GMbh, then we can constrain both rc 
and rc using our observed velocities. From conservation 
of angular momentum for a gas parcel of mass m, we have 


IGMbh 


= m X rcVtrana. 


Thus the BAL gas transverse velocity is Vtrans = 
y/ Ghdrc/r"^ (ignoring any transverse component of its 
velocity away from the black hole across our line of sight 
to the continuum source). The final radial velocity is Vao = 
F^/GMJtl where the scaling factor E’isl.5<F<3.5 
if the wind is accelerated by radiation pressure on ions 
in dust-free gas (see Murray et al. 1995; Laor & Brandt 
2002; Baskin et al. 2014). To solve for rr and rc, we take 
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F = 2.5 ± 1.0 and assume that the observed radial veloc¬ 
ity Vrad,obs equals the terminal velocity Voo- If the latter 
assumption is incorrect, the true vl will be smaller, so we 
call the value we obtain with that assumption rL,max- 
Given the minimum velocity determined for trough A 
above, v > 10,000 kms“^,we find rL,max = TSIjq Rsch = 
0 . 02 t° °2 pc and rc < 186 ± 75 Rsch = 0.04 ± 0.02 pc, 
where the uncertainties on the values of the radii corre¬ 
spond to the values assumed for F. The minimum ve¬ 
locity determined for trough B above was 8,000 kms“^, 
which yields rL,max = 175li®5 Rsch = 0.041q;q 2 pc and 
rc ^ 350 ± 75 Rsch = i3 pc. 


4-1-2 Acceleration 

The above estimate of rL,max assumes that the gas pro¬ 
ducing trough B has reached maximum velocity, which 
may or may not be correct. Here we explore some impli¬ 
cations if that assumption is incorrect. At the small radii 
inferred above, the gas may still be accelerating. We can 
make an order of magnitude estimate of the expected ac¬ 
celeration using some simple assumptions. We stress that 
these assumptions are not unique, only illustrative. 

The radial velocity of a radiatively-accelerated wind 
is approximately v{r) = Uoo(1 — (Murray & Chi- 

ang 1997). The acceleration of the wind is 


a(r) 


dv dv 



If we assume a terminal velocity of Voo = 
60,000 kms“^ for trough B, because the observed veloc¬ 
ity of trough A shows that Civ absorption can be seen 
to that high a velocity, then rc = 78^50 Rsch- If we set 
trough B’s observed velocity Vrad,obs = 40,000 kms“^ = 
v(rc), we find rc = 3.4ri:, = 2651 ^®q Rsch- (Inciden¬ 
tally, that yields a transverse velocity for trough B of 
vtrans = 7, 2 OOI 2100 kms“^, cousisteut with the lower 
limit on the transverse motion velocity we determined 
for a flow-tube in § 4.1). The expected acceleration at 
rc = 265 Rsch for a wind launched at rr = 78 Rsch is 
35111 kms ^ day ^ (the maximum acceleration in that 
model is 86 kms“^ day“^). 

This value is much larger than previous measure¬ 
ments of accelerating BAL winds. For example. Hall 
et al. (2007) measured an acceleration in a C iv trough 
found in SDSS J024221.87-1-004912.6 at approximately 
0.1 kms“^ day“^. The acceleration in J0230, if confirmed, 
would be the largest ever detected in a BAL outflow. 

Using the Gemini South telescope, we have obtained 
a new spectral epoch of J0230 roughly 100 rest-frame days 
after the GEMS epoch of this work. If the above trans¬ 
verse motion variability model is correct, then we predict 
trough B’s centroid velocity will have increased in velocity 
by 3, 500^2400 kms“^ in that data. The results of the new 
observations will be presented in a future paper (Rogerson 
et ah, in preparation). 

This analysis was not done for trough A because we 
have no reliable terminal velocity to suggest the cloud 
might accelerate to. 


4.2 Pure Ionization Parameter Variability 
Model 

In this model, we assume the absorbing clouds are not 
moving across the emitting region of the quasar, and 
thus any variability observed in troughs A and B is due 
to changes in the ionization parameter of the absorb¬ 
ing clouds. In Filiz Ak et al. (2012, 2013), the authors 
observed coordinated variability of distinct Giv BAL 
troughs in the same quasar, even if the troughs are sepa¬ 
rated by as much as 10,000—20,000 kms“^. Other studies, 
such as Grier et al. (2015), observed BAL troughs to vary 
across the entire trough, rather than distinct sections. We 
do not observe either of these behaviours in J0230: we find 
no significant evidence for coordinated variations between 
troughs A and B (they are separated by ~15,000 kms“^), 
and we observe distinct regions of the absorption profiles 
to vary, while others do not (specifically in trough B, see 
§ 2.4.1). Nevertheless, if we assume the changes observed 
in the troughs are due to an ionization state change, we 
can place constraints on the physical properties of the ab¬ 
sorbing gas. Note that in this model only fully saturated 
troughs will not vary. 

The two absorbers responsible for troughs A and B 
cannot have the same distances and densities (includ¬ 
ing density as a function of velocity) to explain the two 
trough’s different responses to the same underlying ion¬ 
izing flux. The exception would be if the absorber closer 
to the quasar significantly reduces the ionizing flux reach¬ 
ing the absorber farther away. Whether the effect is sig¬ 
nificant or not depends on the optical depth to ionizing 
radiation of the absorber closer to the quasar. 

Below, we assume that faster-responding gas has 
higher density. If the changes in trough A are due purely 
to ionization parameter variability, then the high-velocity 
part of this trough has higher density (it responded more 
quickly, and then vanished). If the changes in trough B 
are due purely to ionization parameter variability, then 
the low-velocity part of trough B has higher density (it 
responds faster to ionizing flux changes), and the density 
drops off with increasing velocity. 

One possible pure ionization variability scenario is 
the following. Prior to SDSS, the ionizing flux Fion was 
high, leading to weak absorption. Between SDSS and 
BOSSl, Fion decreased, leading to an increase in Giv ab¬ 
sorption (dense trough A appears). After BOSS2, Fion 
recovered somewhat, leading to weaker trough A absorp¬ 
tion. Between BOSS2 and GEMl, lower-density trough B 
appears in response to the earlier decrease in Fion- The 
above scenario suggests that, barring any major future 
ionizing flux variability, both trough A and trough B will 
decrease in strength with time. Any other trough that ap¬ 
pears will show slower evolution in its EW than trough B 
does, due to the new trough’s required lower density. 


4-2.1 Ionization eonstraints on electron density and 
distanee 

Gonstraints can be placed on the distance from the con¬ 
tinuum source to the absorbing gas, as well as the den¬ 
sity of that gas using the timescale of the variability in 
the absorption. This approach has been used in multiple 
works (see Hamann et al. 1995, Hamann et al. (1997), 
Narayanan et al. (2004), Arav et al. 2012, and references 
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therein). Below, we reproduce the approach taken in Grier 
et al. (2015). 

Consider gas initially in photoionization equilibrium 
in the case where the ionization rate out of ionization 
stage i changes from its equilibrium value li to (1 + /)7i, 
and the rate out of stage i — 1 changes from li-i to 
(1 + /)7i_i,® where / is the fractional change in U. Im¬ 
mediately after this change: 


dm 

dt 


-fmli + fm-ili-i 


-b [— ni{Ii Ri—i) -t- m—iR—i + ni+iRi 


0] (9) 


where the quantity in brackets is the equilibrium value of 
and is therefore zero. In equilibrium, m+ijm = lijRi 
where Ri is the recombination rate to stage i, because ap¬ 
pearance/increase of stage i by recombination from stage 
i + 1 must be balanced by appearance/increase of stage 
i -b 1 by ionization from stage i. Thus we can substitute 
ni_i7i_i = mRi-i = mai-me (using Rmi = ai-iUe, 
where Oi-i is the recombination coefficient to stage * — 1) 
and rewrite as 

^ = - fmli + fmcti-im (10) 

which can be written as 

— = ^ with t* = [-/ (7i - neOi-i)]"^ (11) 

Ui t; 


which is an equation for variations on a characteristic 
timescale t*: mit) = ni(0) exp(t/t*). 

To summarize, for gas which is initially in photoion¬ 
ization equilibrium, the characteristic timescale for den¬ 
sity changes in ionization stage i of some element in re¬ 
sponse to an ionizing flux change can be written as t* 
above (a modified version of Eq. 10 of Arav et al. 2012), 
where — 1 < / < -bc)0 is the fractional change in h, the 
ionization rate per ion of stage i [Ii{t > 0) = (l-b/)7i(t = 
0)], Oi-i is the recombination coefficient to ionization 
stage i — 1 of the ion, and a negative timescale represents 
a decrease in m with time. Note that this equation only 
considers photoionization processes; collisional processes 
are neglected. Gas which shows varying ionic column den¬ 
sities is not in a steady state by definition, but such gas 
can still be in equilibrium with a varying ionizing flux 
if its t* is considerably shorter than the flux variability 
timescale (§ 6 of Pietrini & Krolik 1995). For optically 
thin gas at distance r from a quasar with luminosity Li, 
at frequency v, the ionization rate per ion of stage i is 
given by 


^ /■“ (LRhv)a, 

' A, 4^r2 


( 12 ) 


where is the ionization cross-section for photons of 
energy hv. 

If the absorbing gas is far enough from the quasar 
that li <C rieCy-i-i, then the relevant timescale is Gee = 
IffrieCXi-i (which is just the recombination time of the 
ion in the / = — 1 case where the ionizing flux drops to 
zero), and the observed absorption variability timescale 
constrains the density of the absorber. However, if the 
absorbing gas is close enough to the quasar that R 


® Where we have assumed the fractional change for R and R—i 
is the same. 


rieCti-i, then the relevant timescale is G = —IffR and 
the absorption variations of the ion reflect the ionizing 
flux variations of the quasar, with no density constraint 
derivable just from absorption variations.® An observed 
timescale for variations in optically thin absorption there¬ 
fore constrains the absorbing gas to either have a density 
tie > rimin and r > r^-quai, where requal is the distance at 
which R = nminQi-i, or to be located at r < requal with 
almost no constraint on the density. 

As noted in Arav et al. 2012, there are limitations to 
using timescale arguments to infer physical characteristics 
of an absorber. In that work, the authors indicate ”a more 
physically motivated approach is to use lightcurve simu¬ 
lations that are anchored in our knowledge of the power 
spectrum behaviour of observed AGN lightcurves;” how¬ 
ever, such detailed work is not justified by the relatively 
scarce data available for J0230. 

To determine the constraints on the emergence of 
troughs A and B, we assume a temperature of log T = 4.3 
(Krolik 1999) so that the recombination coefficient is 
aciii = 2.45 X 10“^^ cm® s“^ (from the CHIANTI on¬ 
line database; see Dere et al. 1997, Landi et al. 2013). 
For the simple case of the ionizing flux dropping to zero, 
/ = — 1 and the timescale f* can be approximated as the 
recombination time. Gee ~ XfueOLc iii- 

Using the time between observations of SDSS and 
BOSSl for trough A of 652 days as an upper limit to the 
recombination time, we calculate a lower limit on the den¬ 
sity of the gas to be ne,A 724 cm“®. Using the lower 
limit density of Ue.A -S 724 cm“®, we calculate the mini¬ 
mum distance from the quasar at which that lower limit is 
valid. From its observed flux density at rest-frame 3000 A, 
our quasar has Lhoi — 8.45 x 10^® erg s“^. We adopt the 
spectral energy distribution of Dunn et al. (2010) to cal¬ 
culate Rv Therefore, if the emergence of trough A is due 
to ionization variability, the absorber either has a density 
of rie.A ^ 724 cm“® and is at requal ,a ^ 2.00 kpc, or is at 
r < 2.00 kpc with no constraint on the density. 

Trough B emerged between BOSS2 and GEMl; a 
period of 307 days. Using this as an upper limit to the 
recombination time, we perform the same calculation 
and determine if the appearance of trough B is due to 
ionization variability, the absorber either has a density 
ne,B ^ 1540 cm“® and is at requai,B ^ 1-37 kpc, or is 
r < 1.37 kpc with no constraint on the density. 

Our values for rie are one or two orders of mag¬ 
nitude lower than those found in Grier et al. (2015) 
and Capellupo et al. (2013) (which found values ~ 10® 
cm“® and our values of requal are 10 times larger those 
works (which found values ~100 pc). Further, our values 
of requal are much higher than the launching radius of 
BAL winds expected from theoretical work, which pre¬ 
dict ~ 10“® pc (e.g., Murray et al. 1995). Nonetheless, 
other works have reported outflow radii on similar scales 
to that we infer for J0230 in a pure ionization variability 


® No constraint on rie is derivable even though we can write 
the timescale as 


-faiUe 


( m+i 
V ni 


a»-l Y ^ 
Oi /. 


(13) 


(see equation (2) in Hamann et al. 1997, equation (3) in Arav 
et al. 2015) because in our case rij+i/rij = ncv/^Civ 5 that 
ratio increases more rapidly than rie decreases as the ionization 
parameter increases (Kallman & McCray 1982). 
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model (see Table 10 of Dunn et al. 2010, and references 
therein), and the radius at which a BAL wind is observed 
is not necessarily the radius at which the wind is launched 
(e.g., Faucher-Giguere et al. 2012). 

Finally, we can place an upper limit constraint on 
rie by searching for absorption features from other ions 
of carbon, specifically C ll A1335 A. Given the minimum 
density and r^quai distances determined above, troughs A 
and B are created by absorbers with an ionization param¬ 
eter of Uh — 0.06.^ If we lowered the ionization parame¬ 
ter by a factor of ~50, either by gas at larger radii or at 
higher density, the resulting ionization state would yield 
GII absorption roughly half as deep as the observed G iv 
(see Figure 3 of Hamann et al. 1995). 

A reduction by a factor of ~50 in ionization param¬ 
eter gives us upper limits to both the minimum density 
and the r^q^ai', Therefore, the absorber that caused the 
emergence of trough A has 724 cm~® ^ ne,A ^ 3.62 x 
10"‘(re5„ai,yi/)")^ cm“® and is between r^quai,A ^ r ^ 
7req^ai,A- Similarly, the absorber that caused the emer¬ 
gence of trough B is constrained by 1540 cm“® sC ne,B ^ 
7.70x10 V equal, B IA Cm andrequa/,B ^ I’ ^ ^'l'equal,B ■ 
These upper limits only work in the scenario where we ap¬ 
proximate the recombination time as tree ~ IjneO-c iii- 

In Figure 14, we have plotted the possible values of 
the density of the absorbing gas Ue and the distance the 
absorber is from the source r, given constraints imposed 
by the timescale arguments above for trough B. The ver¬ 
tical and horizontal dashed lines are the locations of the 
requai,B and the minimum electron density rie.B, respec¬ 
tively. Any combination of parameters above the red line 
would have too high a density or too far a distance to 
be ionized to G iv (and lead to the upper limit arguments 
above). There is also a region of too low density or too far 
away that requires too long a timescale for the proper re¬ 
sponse. The allowed regions 1 (between the red and blue 
curves) and 2 (to the left of the green curve) represent 
the combinations of parameters possible. There is also a 
region of too high ionization at low densities and small 
radii which is not visible at the scale shown. Note that 
in the discussion at the end of § 4.2 we assumed that 
faster-responding gas has higher density, although from 
Figure 14 that is only certain if the gas is at r > requal- 
A corresponding plot for trough A would look similar. 


5 SUMMARY 

We have presented the discovery and analysis of two 
extremely high-velocity and highly-variable C iv BAL 
troughs in the quasar SDSS J023011.28-1-005913.6. We re¬ 
trieved 4 spectra of J0230 from the SDSS-I-BOSS archives, 
and obtained 3 of our own spectra using the Gemini 
Observatory. The longest time between observations was 
~ 650 rest-frame days, and the shortest was ~ 10 rest- 
frame days. 

(i) We discovered a Giv BAL trough outflowing from 
J0230 at ~60,000 kms“^ (trough A), the largest veloc¬ 
ity of a BAL wind observed to-date. During follow up 


^ This is for our assumed SED from Dunn et al. 2010, Uh = 
QH At^iahc, with Qh = 6.08 x 10®® hydrogen-ionizing photons 
s“^ and riH = 0.82ne. 



Figure 14. The possible combinations of density and distance 
for the gas that created trough B. The horizontal dashed line 
represents the density if the ionizing flux dropped to zero and 
we use the time between observations as the recombination 
time. The vertical dashed line represents minimum distance 
from the quasar at which the lower limit to the density is valid. 


observations, we discovered a second G iv BAL trough 
outflowing at ~40,000 kms“^ (trough B). See Figure 5. 

(ii) In troughs A and B we observed variability of both 
the depth and shape of the troughs on scales as short as 
10 days in the rest-frame. See Table 2. 

(iii) A dataset of six spectral epochs straddling the 
emergence of both troughs allowed us to rule out some 
simple models of bulk motion as the origin of the variabil¬ 
ity. It also allowed us to propose and test more complex 
and realistic models of bulk motion, such as flow-tube 
geometries and an augmented ‘crossing disks’-|-SS73 sce¬ 
nario. See § 4.1. 

(iv) We found the variability of trough A is best ex¬ 
plained by a ‘crossing disk’ traversing a SS73 disk at 
velocities between 10,000 < v (kms“^) < 18,000. This 
model and velocity range allow us to predict trough A 
will disappear between 350 < t (days) < 1,500 after our 
last observation (GEM3). See § 4.1. 

(v) Trough B is best explained by a flow-tube that has 
recently moved into the line of sight, travelling in the 
velocity range 8,000 < v (kms“^) < 56,000. In this sce¬ 
nario, we have no constraint on how far a flow-tube ex¬ 
tends, and thus cannot predict when trough B will disap¬ 
pear. See § 4.1. 

(vi) Given some simple, conservative assumptions in 
a transverse velocity model, we constrained the distance 
from the black hole to the absorbing gas responsible for 
trough A rc ^ 186 ± 75 Rseh = 0.04 ± 0.02 pc given 
vtrans > 10, 000 kms“^ and for trough B we contrain the 
distance to be rc SS 350 ± 140 Rsch = 0.07 ± 0.03 pc for 
Vtrans > 8,000 hms”*^. See § 4.1.1. 

(vii) If we assume changes to the ionization param¬ 
eter is the reason for the variability observed, the ab¬ 
sorber responsible for trough A either has 724 cm“® ^ 
ne,A ^ 3.62 X W'^{requai,A/r)^ cm“® and is between 
requal,A ^ T ^ 7requal,A, Or is at r < 2.00 kpc with no 
constraint on the density. Similarly, the absorber that 
caused the emergence of trough B is either constrained 
by 1540 cm“® ^ rie.s ^ 7.70 x lV(requai,B/rY cm“® and 
V equal,B ^ r TVequal,B^ Or is at equal ^ 1.37 kpc, or 
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is at r < 1.37 kpc with no constraint on the density. See 
§ 4.2.1. 

Given the resnlts above, we cannot rnle ont bnlk mo¬ 
tion or ionization changes as models of BAL variability. 
More observations of J0230 will, however, allow ns to test 
if our predictions of how the troughs will vary in the fu¬ 
ture are accurate. 
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